Turning a nickelate Fermi surface into a cuprate-like one through heterostructuring 
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Using the local density approximation and its combination with dynamical mean-field theory, we 
show that electronic correlations induce a single-sheet, cuprate-like Fermi surface for hole-doped 1/1 
LaNiOs/LaA103 heterostructures, even though both e g orbitals contribute to it. The Ni 3d 3z 2_ 1 
orbital plays the role of the axial Cu 4s-like orbital in the cuprates. These two results indicate that 
"orbital engineering" by means of heterostructuring should be possible. As we also find strong an- 
tiferromagnetic correlations, the low-energy electronic and spin excitations in nickelate heterostruc- 
tures resemble those of high- temperature cuprate superconductors. 
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The discovery of high-temperature superconductivity 
(HTSC) in hole-doped cuprates [1] initiated the quest for 
finding related transition-metal oxides with comparable 
or even higher transition temperatures. In some systems 
such as ruthenates [2[ and cobaltates [3[ superconductiv- 
ity has been found. However, in these t<i g systems super- 
conductivity is very different from that in cuprates and 
transition temperatures (T c 's) are considerably lower. 

As it became possible to grow transition-metal oxides 
in heterostructures, this quest got a new direction: Novel 
effectively two-dimensional (2D) systems could be engi- 
neered. But which oxides, besides cuprates, are most 
promising for getting high T c 's? 

The basic band structure of the hole-doped cuprates 
is that of a single 2D Cu 3d x 2_ y 2-\ike band which is less 
than half-filled (configuration d 9 ~ h ). In this situation, 
antiferromagnetic fluctuations prevail and are often be- 
lieved to mediate the superconductivity. The Fermi sur- 
face (FS) from this x 2 — y 2 band has been observed in 
many overdoped cuprates and found to agree with the 
predictions of density-functional (LDA) band theory. 

Recently the following idea for arriving at a cuprate- 
like situation in nickelates was presented [4]: Bulk 
LaNiOs (<i 7 ) has one electron in two degenerate e g bands, 
but sandwiching a LaNiOs layer between layers of an in- 
sulating oxide such as LaAlOs will confine the 3z 2 — 1 
orbital in the z-direction and may remove this band from 
the Fermi level, thus leaving the electron in the x 2 — y 2 
band. The possibility of finding bulk nickelates with an 
electronic structure analogous to that of cuprates was 
discarded a while ago [5[ , but heterostructures offer new 
perspectives. 

Indeed, a major reconstruction of orbital states at ox- 
ide interfaces may recently have been observed [6], and 
this kind of phenomenon could lead to novel phases not 
present in the bulk. Extensive theoretical studies of 
mechanisms for orbital selection in correlated systems [7J 
have revealed the complexity of this problem, where de- 
tails of the electronic structure and lattice distortions 
play decisive roles. It is therefore crucial to examine 



nickelate heterostructures by means of state-of-the-art 
theoretical methods and find the optimal conditions for 
x 2 — y 2 orbital selection. 

In this Letter we present results of electronic-structure 
calculations using the merger [8] of density-functional 
(LDA) band theory, which provides an ab initio descrip- 
tion of the materials chemistry, and the dynamical mean- 
field theory (DMFT) [9] , which includes electronic corre- 
lations. We find that the hopping between the x 2 — y 2 
and 3z 2 — 1 orbitals substantially reduces the effects of 
correlations in the 3z 2 — 1 orbital. In this respect, e g elec- 
trons behave very differently than the t^g electrons, which 
have no inter-orbital hopping on a square lattice. Never- 
theless, we do find that the correlations may sufficiently 
shift the bottom of the hybridizing e g bands relatively to 
each other to yield a FS with only one sheet. This sheet 
has predominantly x 2 — y 2 character and a shape like in 
the cuprates with the highest T c max (T c at optimum hole 
doping) [10[, but even more extreme. Moreover, stretch- 
ing the in-plane lattice constants by suitable choice of 
substrate reduces the correlation-strength needed to pro- 
duce a single-sheet FS. Since we also find strong anti- 
ferromagnetic fluctuations, somewhat larger than in the 
cuprates, nickelate heterostructures hold the basic ingre- 
dients for high-temperature superconductivity. 

Here we give results for the simplest, 1/1 superlattice 
LaNi0 3 /LaA10 3 = LaO-Ni02-LaO-A102 shown in the 
left-hand side of Fig. 1. For the in-plane lattice constant 
a we first took that of SrTi03, often used as substrate, 
whereby the Ni-0 and Al-0 distance in the x- and y- 
directions became: £Ni-o = 1-95 A, not far from the value 
in pseudo-cubic LaNiOs. The lattice constant c we took 
as the sum of those of pseudo-cubic LaNiOs and LaA103, 
whereafter the position of apical O was relaxed within 
the LDA [11] to yield: z m . = 1.91 A, i.e. 2% smaller 
than XNi-o- Next, we expanded the LaNiOa/LaAlOs het- 
erostructure in the x-and ^/-directions by 3%, as might 
be achieved by growing LaNiOs /LaAlOs on a PrScOs 
substrate, to yield x^-o = 2.01 A. With the concomi- 
tant 6% contraction in the z-direction, relaxation of the 




FIG. 1: The 1/1 LaNi0 3 /LaA10 3 heterostructure (left) and its LDA (iVMTO) bandstructures without (center) and with (right) 
strain. The Bloch vector is along the lines T (0, 0, 0) - Z (0, 0, f ) - R (0, f , f ) - A (f , f , £ ) - Z (0, 0, f ) . The shading gives 
the x 2 — y 2 vs 3z 2 — 1 e g Wannier-function character. 



apical-oxygen position within the LDA finally lead to: 
2 N i_o = 1.81A. 

Fig. 1 shows the LDA energy bands for the two differ- 
ently strained heterostructures in a 5-eV region around 
the d 7 Fermi level (=0) . The two solid bands are the 1/4- 
full Ni-0 pda antibonding e g bands, which are pushed 
up above the less antibonding Ni-0 pdir t^g bands (thin 
bands) lying below —1 eV and well above the Ni-O, Al-O, 
and La-0 bonding bands below the frame of the figure. 
The antibonding Al-0 and La-0 bands (thin bands above 
1-2 eV) lie respectively 9 and ~ 5 eV above their bond- 
ing counterparts, and as a result there is a comfortable 
2-3 eV gap above the top of the antibonding t^ g bands 
in which the two antibonding e g bands reside. 

The shading (coloring) of the e g bands gives the rela- 
tive x 2 —y 2 and 3z 2 — l characters in the Wannier-function 
representation of these two bands, as calculated with the 
iVth-order muffin-tin-orbital (iVMTO) method and N=2 
[l2[. We see that in the "nodal" k x =k y plane containing 
the TZ and AZ lines the x 2 - y 2 (\m\ =2) and 3z 2 - 1 
(m=0) Wannier orbitals cannot mix [13[. The bottoms 
of both bands are along TZ, i.e. for k x =k y =0. That of 
the x 2 — y 2 band is at —1.5 eV and does not disperse with 



—3cosak independently of k z because tdd5 is negligible. 
This means that both bands extend from — 3 \tdda\ at 
(0,0,0) to +3 \t d da\ at (I , f , f ) in the bulk. Substitut- 
ing now every second LaNiOs layer by an "insulating" 
LaA103 layer, forces the Bloch waves to have nodes ap- 
proximately at the AIO2 planes, so that only waves with 
kz ~ J~ ~ ~ are allowed. As a consequence, the bottom 
of the 3z 2 — l band is pushed up by ~2 \tdda\ > i- e - by ~l/3 
the e g bandwidth. The exact position of the nodes, and 
hence the upwards shift of the 3z 2 — 1 band, depends 
on the scattering properties of the insulating layer. This 
suggests that the band structure can be tuned by choice 
of the insulating layer. 

A further factor influencing the orbital separation is 
the tetragonal Jahn- Teller (JT) distortion of the nickel- 
centered oxygen octahedron. Since the x 2 — y 2 and 3z 2 — l 
Wannier orbitals antibond with oxygen, flattening the oc- 
tahedron (zNi-o < #Ni-o) moves the energy of the former 
orbital down, and that of the latter up. However, this 
crystal-field splitting is little effective in achieving orbital 
separation for configuration d 7 ~ h because the e g Bloch 
sums at the bottom of the cubic band at (0,0,0) have 
no oxygen character, so only energies higher up in the 



k z , while that of the 3z 2 — 1 band is at —0.5 eV at T and e g band are effected. For JT-nattening to be effective, 



disperses upwards to —0.1 eV at Z. The bottom of the 
3z 2 — 1 band is thus 1 eV ^ 1/4 e 5 bandwidth above 
that of the x 2 — y 2 band. Straining by 3% is seen to shift 
the bottom of the 3z 2 — 1 band up by further 0.2 eV. 
The LDA FS thus has two sheets, and reducing it to one 
would require moving the 3z 2 — 1 band above the x 2 — y 2 
band at T by an additional 0.5 eV for the unstrained and 
by an additional 0.3 eV for the strained super lattice. 

That the x 2 — y 2 Wannier orbital is more populated 
than 3z 2 — 1 (the ratio is 70/30 for the unstrained su- 
perlattice) is mainly due to the confinement in the z- 
direction. Consider for simplicity the dispersions in the 
k x =± k y =k planes where the 3z 2 — 1 and x 2 — y 2 orbitals 
do not hybridize: In cubic, bulk LaNiOs, e^^-i (k,k z ) « 
— cos ak — 2 cos ak z , with respect to the center of the 



e g band and in units of \tddo 
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confinement is therefore a prerequisite. This is clearly 
seen from the LDA bands for the 3% strained superlat- 
tice on the right-hand side of Fig. 1: Whereas the strain 
moves the top of the x 2 — y 2 band down and that of the 
3z 2 — 1 band up, the bottom of the x 2 — y 2 band is not 
affected, and that of the 3z 2 — 1 band is shifted up only 
because it has antibonding oxygen character correspond- 
ing to k z ~-^ rather than to k z =0. 

For the undoped (<i 9 ) cuprates, the LDA bandstruc- 
tures are roughly similar to this, but the antibonding 
3z 2 — 1 band is now full and lies in the region of the t^ g 
bands. Filling this band has annihilated the pda bond 
to apical oxygen and thereby caused zqu-o to increase 
well beyond xc u -0, whereby the antibonding push-up of 
the 3z 2 — 1 band has been lost. The half- full pda an- 
tibonding x 2 — y 2 band lies a bit lower with respect to 



the O and cation bands than in the nickelates because 
the position of the 3<i-level in Cu is lower than in Ni. 
However, the shape of this cuprate conduction band near 
half filling is not unlike that of the lowest e g band in the 
nickelate heterostructures, in particular for the cuprates 
with the highest T c max . Specifically, LDA calculations 
for a large number of cuprate families have revealed that 
whereas the dispersion along the nodal direction (ZA) is 
always the same, the energy of the saddlepoints at (^, 0) 
and (0, -) , i.e. at R, depends on the material and is 
lower for materials with higher T cmax [10]. The reason 
for this correlation is not understood, but the reason for 
the change of band shape is clearly that the x 2 — y 2 orbital 
is hybridizing with a material-dependent axial (\m\ =0) 
orbital whose energy lies ~10 eV above the Fermi level, 
but falls for cuprates with increasing T c max . This axial 
orbital is essentially the antibonding linear combination 
of Cu 4s and apical O 2p Zl so that its energy falls if their 
interaction decreases, e.g. by increasing zc u -0- Concomi- 
tant with this change of band shape is a concentration 
of the conduction-band Wannier function onto the Cu02 
layer, away from the perpendicular direction. Instead of 
using the energy of the axial orbital as band-shape pa- 
rameter, one uses a dimensionless parameter, r, which 
for materials with low T c max (< 50K) becomes the ratio 
t'/t of the 2nd to the lst-nearest-neighbor hopping in- 
tegral. The cuprates with the highest T c max (~ 140 K) 
have r ~ 0.4. If one could lower the energy of the axial 
orbital right down to the Fermi level, r would have the 
value 1/2. 



band Hubbard Hamiltonian: 
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Here, the on-site Coulomb terms, namely the intra and 
inter-orbital Coulomb repulsions, U and V = U — 2 J, 
as well as the Hund's exchange, J, have been added 
to the LDA e g Wannier- function Hamiltonian, H^ m ,. 
This Hubbard Hamiltonian we solve for 1/4 filling in the 
single-site DMFT approximation for the paramagnetic 
phase and at a temperature so high (1160 K = 0.1 eV/k B ) 
that we can afford using the Hirsch-Fye Quantum Monte 
Carlo method. 
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FIG. 2: Cross-section of the LDA (left) and LDA+DMFT 
(U — 6.7 eV) FS with the k z —0 plane for the unstrained 1/1 
heterostructure. e^-orbital characters are coded as in Fig. 1. 



This axial-orbital model also applies to the e g bands of 
nickelate heterostructures, but due to the short distance 
to apical oxygen, the axial orbital is now essentially the 
antibonding linear combination of Ni 3d 3z 2_ 1 and apical 
O 2p z . Its energy is that of the 3z 2 — 1 band at T, and 
since this is below ef for the LDA bands shown in Fig. 
1, they have r > 1/2. Engineering these heterostructures 
should presumably first aim at reducing r towards that 
(~0.4) of the cuprates with the highest T c max , i.e. at 
moving the energy of the second band at T well above 
Ef- This requires increasing the interaction between Ni 
3d% z 2 _i and apical O 2p Zl e.g. by reducing ZNi-o- As we 
shall see, this is helped by the electronic correlations, but 
does not necessarily lead to HTSC, because although the 
same value of r gives the same band shape for nickelates 
and cuprates, their conduction-band Wannier orbit als are 
not identical. 



Having studied the materials dependence of the LDA 
band structures, and having found that the conduction 
bands in the paramagnetic phase are well separated from 
all other bands, we can study the effects of Coulomb cor- 
relations in the nickelate heterostructures using the two- 



Our DMFT calculations confirm the common expec- 
tation that, for a metallic multiband system, the main 
effect of the Coulomb correlations is to enhance the split- 
ting between the subbands such as to reduce the density 
of states at the Fermi level. Specifically, for the undoped 
superlattice with J=0.7eV and U increasing, we find that 
the bottom of the 3z 2 — 1 band is driven up and passes 
the Fermi level when U exceeds 6.4 eV for the unstrained 
and 5.7 eV for the strained structure. Hereafter the FS 
has only one sheet, a large (^, ^-centered hole cylinder 
whose shape can be seen from Fig. 2 to be similar to that 



found in the cuprates with the highest T c 



but even 



more extreme. It is of course possible that the strong 
nesting of this FS makes it unstable with respect to spin 
or /and charge-density waves with q x ~^ an d Qy^j^i smi_ 
ilar to what has been found in cuprates. At the point 
where the second sheet disappears, r=l/2 and the ratio 
between the x 2 — y 2 and 3z 2 — l populations has increased 
to 80/20 for the unstrained - and beyond for the strained 
- superlattice. Reasonable changes of J slightly influence 
details of the Hubbard subbands, but not the physics of 
the transition. 

The remaining half-full band undergoes a Mott- 
transition when U exceeds 7.4 eV for the unstrained and 
6.5 eV for the strained superlattice. For comparison, a 



half-full cuprate band undergoes a Mott transition in 
DMFT for a critical value of U which increases with r 
and takes the value 4.5 eV for r=0.4 [14]. This behav- 
ior for the cuprates is thus in line with what we find for 
the nickelate heterostructures where tlda (unstrained) > 
^lda (strained) ~ 1/2, and this supports our hope that 
the nickelates can be engineered such that, like in the 
cuprates, hole-doping will suppress the Mott transition 
and produce superconductivity. For nickelates there is 
even the possibility of engineering the e g bands such that 
the real value of U falls between the one needed to reduce 
the FS to a single sheet and the one needed to eliminate 
this sheet by a Mott transition. If this can be achieved, 
superconductivity in the nickelates may occur even with- 
out doping. This is a remarkable result. 
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FIG. 3: Energy levels for the unstrained two-site model with 
U — 6.4 eV as a function of the splitting A between the 
energies of the 3z 2 — 1 and x 2 — y 2 Wannier orbitals. The 
LDA value of A is indicated by the dashed line. Of (Oaf) 
denotes a configuration with the same (different) orbital(s) on 
the two sites. 



Next, we need to estimate the strength of antifer- 
romagnetic correlations, which are believed to play a 
central role in the physics of the cuprates. Since our 
LDA+DMFT calculations would be prohibitively expen- 
sive for the study of low temperature magnetic proper- 
ties, we merely diagonalized the two-site version of the 
Hubbard Hamiltonian obtained by Fourier transforma- 
tion of i^ m / and truncation to a diatomic molecule di- 
rected along x. The energy levels are presented in Fig. [3] 
as functions of the the difference, A, between the ener- 
gies of the 3z 2 — 1 and x 2 — y 2 Wannier orbitals. The 
ground state is always a spin singlet. Increasing A from 
to oo leads to demixing such that the orbital configura- 
tion changes from 3x 2 — 1 to x 2 — y 2 . For the LDA value, 
the orbital character is already close to x 2 — y 2 . From the 



distance between the singlet ground state and the triplet 
first excited state, we estimate the magnitude of the anti- 
ferromagnetic coupling constant to be J af ~ 0.2 eV, i.e. 
somewhat higher than in cuprates. 

Altogether, our analysis of the 1/1 LaNiOs/LaAlOs 
system shows that heterostructuring of d 7 nickelates is 
promising because their physics contains the main ingre- 
dients of high-temperature superconductivity. In partic- 
ular, we find that electronic correlations reduce the FS 
to a single sheet whose shape is similar to the one in 
the hole-doped cuprates with the highest T cmax . This 
sheet has not only x 2 — y 2 , but also 3z 2 — r 2 charac- 
ter, and this gives a new twist to the intensive discussion 
of orbital-selective Mott-Hubbard transitions. Substrate- 
induced strain and/or use of insulating layers different 
than LaA103 may tune the FS shape and may enable 
superconductivity without doping. 

Discussions with J. Chakhalian, H.-U. Habermeier, 
and T. Saha-Dasgupta are gratefully ackowledged. 
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